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This review summarizes recent attempts to reconstruct the expansion history of the Uni- 
verse and to probe the nature of dark energy. Reconstruction methods can be broadly 
classified into parametric and non-parametric approaches. It is encouraging that, even 
with the limited observational data currently available, different approaches give con- 
sistent results for the reconstruction of the Hubble parameter II[z) and the effective 
equation of state w{z) of dark energy. Model independent reconstruction using current 
data allows for modest evolution of dark energy density with redshift. However, a cos- 
mological constant (= dark energy with a constant energy density) remains an excellent 
fit to the data. Some pitfalls to be guarded against during cosmological reconstruction 
are summarized and future directions for the model independent reconstruction of dark 
energy are explored. 

1. Introduction 

The accelerated expansion of the universe has now been confirmed by several in- 
dependent observations including those of high redshift type la supernovae, and 
the cosmic microwave background (CMB) combined with the large scale structure 
of the Universe i'2.3^4,5^ Another way of presenting this kinematic property of the 
Universe is to postulate the existence of a new entity - dark energy (DE). This 
latter statement is dynamical in nature and therefore requires some assumptions to 
be made about the form of gravitational field equations governing the evolution of 
the (observed part of the) Universe. 

Although observationally well established, no single theoretical model provides 
an entirely compelling framework within which cosmic acceleration or DE can be 
understood. Indeed, the very many models of DE existing in the literature illustrate 
that its nature is still very much an enigma. At present, all existing observational 
data are in agreement with the simplest possibility of DE being a cosmological 
constant A with pA = A/SttG = const ~ 10~*^GeV'* (inside ^ 2a error bars in 
the worst case).** This case is internally self-consistent and non-contradictory. The 
extreme smallness of the cosmological constant expressed in either Planck, or even 

^/i = c = 1 is used throughout the paper. 
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atomic units means only that its origin is not related to strong, electromagnetic and 
weak interactions (in particular, to the problem of the energy density of their vac- 
uum fluctuations).'^ Although in this case DE reduces to only a single fundamental 
constant we still have no derivation from any underlying quantum field theory for 
its small value. 

Within this context it is but natural that other possibilities admitting a (slightly) 
variable dark energy have also been actively studied by the scientific community in 
recent years. Moreover, it is interesting that properties of the currently observed 
'late DE' are qualitatively similar to those of an 'early DE' which is believed to have 
given rise to accelerated expansion (inflation) in the early Universe. However in the 
case of the latter, there are sufficient grounds to support the view that 'early DE' 
was unstable and, thus, more complicated than a cosmological constant. So, it is 
natural to conjecture by analogy that the same might also be true of 'late DE'. 

DE models proposed to account for the present cosmic acceleration include: 

(i) Quiessence with w = pde/pde = constant, the cosmological constant A (w = 
— 1) is a special member of this class. 

(ii) Quintessence models which are inspired by the simplest class of inflationary 
models of the early Universe and employ a scalar field rolling down a potential 
V{(j)) to achieve late-time acceleration. Quintessence potentials with V"V/ (V)^ > 1 
have the attractive property that dark energy approaches a common evolutionary 

'tracker path' from a wide range of initial conditions. 

(iii) The Chaplygin gas model (CG) has the equation of state p oc — 1/p and evolves 
as p = A-\- B(\ -\- zY" where z is the redshift, z = a{to)/a{t) — 1. It therefore 
behaves like dark matter at early times (z 1) and like the cosmological constant at 
late times. CG appears to be the simplest model attempting to unify DE and non- 
baryonic cold dark matter. 

(iv) 'Phantom' DE (w < -1). 

(v) Oscfllating DE. 

(vi) Models with interactions between DE and dark matter. 

(vii) Scalar- tensor DE models. 

(viii) Modified gravity DE models in which the gravitational Lagrangian is changed 
from R to F{R) where R is the scalar curvature and F is an arbitrary function. 

(ix) Dark energy driven by quantum effects. 

(x) Higher dimensional 'braneworld' models in which acceleration is caused by the 

leakage of gravity into extra dimensions. 

(xi) Holographic dark energy, etc. 

See the reviews ^'"^ for an exhaustive list of models and references. However, none 
of these models leads to the reduction of the number of fundamental constants 
(parameters in the microscopic Lagrangian) as compared to standard ACDM. In 

However, the empirical relation pA ~ where m^y is some characteristic neutrino rest-mass (the 
lightest one ?) suggests that vacuum energy of the intera<;tion responsible for non-zero neutrino 
rest-masses may be relevant for a non-zero A. 
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other words, at the current state-of-the-art, DE requires at least one new parameter 
whose value is set from observations. 

Models with variable DE can be broadly divided into two main classes: 

1 . Physical DE, in these models DE is the energy density of some new, very 
weakly interacting physical field. 

2 . Geometrical DE (otherwise dubbed modified gravity models). In these mod- 
els the gravity equations do not coincide with those of Einstcinian general relativity. 
However, it is usually possible to re- write the new equations in the conventional Ein- 
steinian form by transferring all additional terms from the l.h.s. into the r.h.s. of the 
Einstein equations and referring to them as an effective energy-momentum tensor of 
DE (see Sec. 2 below). (Precisely this happened to the cosmological constant which 
originally appeared in the l.h.s. of Einstein's field equations but is now felt by many 
to constitute an effective matter term such as vacuum energy, etc.) 

Another category, which is even more important from the observational point 
of view, arises in response to the question whether or not the description of DE 
requires a new field degree of freedom (= a new kind of matter). If the answer is 
in the affirmative then DE may be considered as being 'induced' by other kinds of 
matter. All physical DE models and many geometrical ones belong to this category 
but there do exist geometrical DE models which do not (for instance the F{R) 
model with the Palatini variation of its action). 

Faced with the increasing proliferation of DE models each with its own physical 
motivations and assumptions, a concerned cosmologist can proceed in either of two 
ways: 

(i) Test every single model against observations. 

(ii) Try and ascertain properties of dark energy in a model independent manner. 
In this article, we proceed along route (ii) and attempt to review both the 

successes as well as difficulties faced by methods attempting to reconstruct the 
properties of dark energy directly from observations in a model independent manner. 



2. Model independent reconstruction of Dark Energy 

Before attempting to determine its properties, we first need to provide a definition of 
dark energy. A traditional approach is to use the Einstein form of the gravitational 
field equations 

R^. - \9^uR = SttG T^t^ + T^J"^ (1) 

as providing a definition of the effective energy- momentum tensor T^^f of DE. Here, 
the summation over a in the r.h.s. includes all types of matter known from laboratory 
experiments (protons, neutrons, photons, neutrinos, etc.) as well as non-relativistic 
non-baryonic cold dark matter (whose energy-momentum tensor is dust-like in the 
first approximation 0<p<^p).G = Gq = const is the present value of Newton's 
gravitational constant. 
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This definition has a number of advantages: 

(i) it is simple, well-defined and self- consistent; 

(ii) it treats physical and geometrical DE on an equal footing; 

(iii) in the absence of direct physical interaction between DE with known forms of 
matter or with cold dark matter, the DE energy-momentum tensor is conserved: 

One should stress here that using (1) we automatically ascribe terms describing 
(possible) gravitational interactions between DE and non-relativistic matter, as well 
as the matter energy-momentum tensor multiplied by a change in the effective 
gravitational constant, to the DE energy- momentum tensor. The latter possibility 
arises, for instance, in scalar-tensor DE models. 

When applied to a homogeneous and isotropic Friedmann-Robertson- Walker 
(FRW) cosmological model, eqs. (1) reduce to two algebraically independent equa- 
tions: 

jj2 



Pa + PDE^ - , (2) 
a 3 yY.^Pa + ^Pa) + PDE + SpDEj (3) 



3 

a AttG 



where a{t) is the FRW scale factor and the Hubble parameter H{t) = a/a. 
At late times when radiation may be neglected, one gets 

PDE = ^(1-J^J (4) 

where we have omitted the contribution from the curvature term in (2) for simplicity. 

is the total density of non-relativistic matter in terms of its critical value. 
Similarly, the expression for the deceleration parameter 

q = -a/aH'^ = ^„^^} x-1, x = l + z (5) 
H{x) 

(prime implies differentiation with respect to x) takes the form: 

1 

Dividing (6) by (4) we get the following expression for the effective equation of 
state (EOS) of dark energy {iv = Pde/pde) '■ 

2q{x) - 1 _ {2x/i) dlnH / dx-1 
"^^^^ ~ 3{l-n^{x)) = 1 - {Ho/HY^mQ x^ ■ ^ ^ 

For physical DE, the EOS makes physical sense. However, this is not so for geomet- 
rical DE for which the acceleration of the Universe is caused by the fact that the 
field equations describing gravity are not Einsteinian. 
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As a specific example of geometrical DE consider the widely studied DGP 
braneworld ^ for which 



where Ic = m'^ /M^ is a new length scale and m and A-I refer respectively to the four 
and five dimensional Planck mass. The acceleration of the universe in this model 
arises not because of the presence of DE but due to the fact that gravity becomes 
five dimensional on length scales R > /<; = 2Hq^(1 — Om)^^- The contrast between 
(2) and (8) makes it abundantly clear that for models such as DGP the EOS in (7) 
is an effective quantity {w = Weff), which may still be useful for descriptive purposes 
but which no longer represents any fundamental physical property of an accelerating 
universe. Indeed, instances are known when WcS < — 1 even when matter itself sat- 
isfies the weak energy condition p + P > S-iO'H. It may be instructive to note that, 
for geometrical DE, w{z) may show pathological behaviour in certain cases, such 
as the presence of poles at finite values of redshift, w{zp ±oo), even though the 
underlying cosmological model is completely well behaved (sec for instance ^^'^^). 
For such models, the deceleration parameter q{z) and other geometrical parameters 
prove to be more robust for determining DE properties than the EOS. We shall 
return to these important issues in Section 6. 

Observational tests of DE rely on an accurate measurement of at least one of 
the following quantities: 



(1) The luminosity distance 



Dl{z) _ r dz' 
l + z 



(2) The angular size distance 



l + z Jo H{z') 

(3) The coordinate distance 



Da{z) = —-I (10) 



(11) 



r{z) = r 

Jo 



dz' 



H{z') 

Here, we have assumed that the universe is spatially flat for simplicity. Note that 
the distance duality relation Dl = (1 + zYDa{z) which follows from (9) and (10) is 
valid only for metric theories of gravity Its violation (if observed) could, therefore, 
be used to probe alternative theories of gravity. Existing data, however, appear to 
support it (within observational errors) . 

In all of the above expressions the value of the Hubble parameter can be 're- 
constructed' through a relation similar to the one given below for the luminosity 
distance I6,i7,i8,i9,20. 



H{z) 



d_ f DL{z) 
dz\l + z 



-1 



(12) 
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Differentiating a second time allows one to reconstruct the equation of state of 
DE (7). Equations (7) and (12) immediately inform us that iu{z) will be a noisier 

quantity than H(z) since two successive differentiations are needed for the recon- 
struction Dl w{z) while a single differentiation suffices for Dl H{z). This has 
led several people to suggest H{z) (or /9de(-2)) as being better suited for providing 
a model independent description of properties of DE. Another important difference 
between H{z) defined in (12) and w{z) in (7) is that the former is independent of 
the value of the matter density parameter fi^ while the latter is not. As a result, 
uncertainties in the current value of affect the reconstruction of the EOS far 
more profoundly than they do H{z). We shall return to this issue in section 5. 

Knowing H{z) (either through (12) or using corresponding relations for other 
observational tests), allows us to extend cosmological reconstruction to other im- 
portant physical properties of the Universe including: 



Its age 



The deceleration parameter (5) and the equation of state (7); 
the electron-scattering optical depth to a redshift Zreion 

r(^.e.onj {1 + Z)H{Z)' ^'^> 

where Tig is the electron density and ctt is the Thomson cross-section describing 

scattering between electrons and CMB photons. 

The product dA{z)H{z), which plays a key role in the Alcock-Paczynski 
anisotropy test . 

• The product (i^(z)iJ~^(z), which is used in the volume-redshift test 

• The parameter A associated with the determination of the baryon acoustic peak 



23 

2/3 



^ _ vHrnO 



/l(^l)V3 



Zl Jo 



dz 



0.469 ±0.017, (15) 



h{z) 

where h{z) = H{z)/Ho and zi = 0.35 is the redshift at which the acoustic scale 

has been measured in the redshift sample. 
• The 'shift' parameter R associated with the CMB '''24.25 

dz_ 
Mz) 

For quintessence, it is additionally possible to reconstruct its potential, since the 
Einstein equations 

if = -47rG(p„ + </.2) , (17) 



R=^/q^ ^ = 1.7 ±0.03. (16) 
Jo 
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can be rewritten as (see, for instance, i^'^^): 

^J{x) = ^ - - -n^ox , (18) 

87rGfdiy_^dlnH_n^ 

Integrating (19), one determines 4>{z) to within an additive constant. The inversion 
(j){z) z{(p) followed by substitution into (18) allows us to reconstruct V((j)), since 
H{z) and its first derivative can be determined from observations using (12). The 
presence of rimo in (18) and (19) implies that the value of the matter density must 
be known rather precisely for an accurate reconstruction of V{4>). 

In several important cases (18) and (19) have led to a closed form expression for 
V {(/)), for instance: 

• DE with a constant equation of state — 1 < iz; < is described by ^'^^ 

= '^^^''"-l^oa""""^'^'"" sinh-^" I Ha/t^('/'-<^o + <^i) 1 ,(20) 



IBttGO^o 




where 



0. = ——, (f>Q = (f>{to), <^i = W— — — In — = . (21) 

\w\ V dttG \w\ V^^mO 

Consequently, a universe consisting of such a scalar field will have expansion 
properties which closely mimic a different kind of DE, e.g., a tangled network 
of cosmic strings (w; = —1/3) or domain walls [w = —2/3). 
• The Chaplygin gas model which unifies dark matter and DE and has p = —A/p 
can be described by the scalar field potential (in the absence of additional 
cold dark matter and neglecting baryons and photons): 

VW = (co.l,(2v^iG« + ^^^^l^^.) . (22) 

Note that the behaviour of the Chaplygin gas may also be modelled completely 
differently using a scalar field with the Born-Infeld kinetic term ^^'^^ (a specific 
type of k-essence). This illustrates once more that, even for physical DE, the 
equation of state w{z) does not uniquely define an underlying field-theoretical 
model (see also ^^). 

From the fact that the left-hand side of Eq. (19) is always non-negative, follows 
an important restriction on the expansion law for the Universe which must be 
satisfied before attempts to reconstruct the potential are made, namely 

^>m^oH^{l + zf . (23) 
ctz 



■^The Born-Infcld Lagrangian C = — Vby 1 — <t>,ij,(i>'^ is quite different from the Quintessence La- 
grangian C= 5</>^ — V {<!>). 
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Equation (23) is simply a restatement of the weak energy condition + p,p > 0. 
Integrating (23), we get the relation 

H\z)>H^,{l + n^o{l + zf) (24) 

which is easier to verify observationally. 

In the case of physical DE, with its implied minimal coupling to gravity, the 
expansion history H{z) can also be reconstructed from the growth rate of inhomo- 
gencous density perturbations in the non-rclativistic matter component on scales 
significantly less than the Hubble radius H~^, (provided perturbations in the DE 
component can be neglected and the effective gravitational constant does not change 
with time) Indeed, the linearized perturbation equation 

S + 2HS - ATrGpmS = (25) 

can easily be inverted, with the result: 



H{z) = H{0) 



(26) 



Furthermore, an interesting relationship follows between the current value of the 
matter density ilmo and 5{z): 

nmO = ^"(0)(^3^~^rfz^ , (27) 

which could provide a consistency check on direct observational determinations of 
rimo- At sufficiently low redshifts z < 1, the advent of deep redshift surveys probing 
large scale structure may help in reconstructing H{z) using (26). Although it is 
unlikely that the value of 5{z) for z > 1 will be reliably known in the near future, 
the fact that the Universe is expected to become matter dominated at fairly low 
redshifts, ilm I aA, z ^ 1, allows us to use the spatially flat matter dominated 
solution ^ oc {1+ z)~^ to extrapolate to higher redshifts and thereby evaluate (27). 

As a result, for physical DE we now have two independent methods to recon- 
struct the same quantity H{z)/Hq: (i) through quantities like Dl{z) referring to 
an unperturbed FRW background, and (ii) through 6{z) describing the growth of 
perturbations on small scales. It therefore follows that one should be able to recon- 
struct the density perturbation 5(z) from observations of the luminosity distance. 
Consider the quantity E{z) = HqDl{z)/{1 + z) which can also be written as the 
difference of conformal times: E{z) = aoHo{r]{Q) — r]{z)), r] = J dt/a{t). For our 
purpose it will be useful to invert the quantity E{z) (determined from observations) 



Galaxy peculiar velocities produced by these matter inhoinogeneities result in small angular 
dependent corrections to the luminosity distance Dj^(z) and other quantities referring to the 
evolution of a FRW background that, in turn, can be used to determine H{z) 32,33,34^ 
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to z{E). Next, we again assume that DE does not cluster and rewrite (25) in the 
following integral form: 



S{E{z)) = S{0)+6'{0) / {l+z{E))dE+^ / {l+z{Ei))dEi / S{E2)dE2.{28) 



Here 6' = d6/dE = dS/dz since E{z) « 2 for 2 <C 1. Clearly (28) can be solved 
iteratively if the value of 5'{0) is known. The value of (5(0) does not really matter 
since we are only interested in the ratio 6{z)/6{0). The requirement 5{z = 00) = 
leads to the integral condition: 



and, as in the case of (27), one needs to postulate some reasonable interpolating 
behaviour for 6{E{z)) at large redshifts for which no supernova data exist at present. 
Since Eiss = Jq^"" dz/h{z), the value of this quantity should be possible to determine 
from the CMB. 

Both (26) and (28) can be used as consistency checks for physical DE which, 
by assumption, is minimally coupled to gravity. For geometrical DE on the other 
hand, the linearized perturbation equation is usually modified from its conventional 
form (25), see for instance ^ for scalar-tensor gravity or for the DGP model. 
As a result, DE reconstruction becomes more complicated. In particular, a generic 
scalar-tensor model of DE depends on two arbitrary functions V{(j)) and F{(l)) - 
the scalar field potential and its coupling to gravity respectively. In this case one 
needs information from both Dl{z) and 5{z) in order to reconstruct the DE model 
unambiguously ^ . 

To summarize, we have shown in this section that the reconstruction of the 
expansion rate H{z) and other DE properties from observational data is well defined 
and unambiguous from the mathematical point of view. In practice, however, the 
situation is much more diflacult since all observational functions such as Dl{z) 
and S{z) are noisy and known only at discrete values of the rcdshift {zi, Z2....zn} 
(associated with the redshifts of N supernovae in the case of Dl{z)). Thus, it is 
impossible to directly difl[erentiate them with respect to redshift as formulae (12) 
and (26) require. Add to this the fact that the dispersion in the luminosity distance 
is not expected to get significantly better than crindt = 0.07, and one is confronted 
with a noisy quantity DL^Zi) sampled at a set of discrete intervals {zi}. Therefore, 
to convert from D^lzi) to the function H{z) defined at all redshift values within the 
interval (say) < 2; < 2 using (12) clearly requires a crucial additional step: some 
sort of smoothing procedure. This is usually accomplished using either parametric 
or non-parametric reconstruction which we turn to next. 







(29) 



•^Thc fact that in this case 5(z) is expected to be different from the reconstructed expression (28) 
provides an interesting possibihty to distinguish between physical and geometrical models of DE. 
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3. Parametric reconstruction 

This approach is based on the assumption that the quantities Dl{z), H{z),w{z) 
vary 'sufBciently slowly' with redshift and can therefore be approximated by a fitting 
formula (ansatz) which relies on a small number of free parameters ai, i = 1,N. 
The ansatz for DL{z,ai),H{z,ai),w{z,ai) is compared against observations, and 
the values of the free parameters Oj are determined using a minimization procedure 
(usually maximum likclyhood). Quite clearly a successful fitting function should, 
in principle, be able to faithfully reproduce the properties of an entire class of DE 
models. Different parameterizations have been suggested for: Dl 16,17,19,20,37^ H{z) 

26,38,39,40 ^^^^ 41,42,43,44,45,46,47,48,49,51,52,53,54,55 ^^^^ V{z) ^^'^"^ 

We attempt to summarize some of these approaches below. 

(1) Fitting functions to the luminosity distance Dl 

An approximation to a function can easily be generated by expanding it in a 
Taylor series about some redshift zq- Applied to the luminosity distance (with 
zq = 0), this method gives the following fitting function 

^ = f:a.z^. (30) 

Unfortunately, it appears that in order to accurately determine quantities of 
interest such as H{z) and w{z) one must take a large number of terms in (30) 
which greatly increases the errors of reconstruction '^^ and reduces the reliability 
of this ansatz. 

A more versatile Pade-type ansatz was suggested in 



HoDl{z) 

1 + z 



Aly/^- 1+Al 



A2x + AsV^^ + 2-A,-A,^As\ ' ^ + ^^^^ 

which is able to exactly reproduce the results both for CDM (fi^ = 1) and the 
steady state model (f2A = !)• 

Another accurate fit with a greater number of free parameters is 

^^=r?(l) -.(,), (32) 

7^{y) = 2a [y-^ + (iy-^^ + jy-* + 5y-^ + a] "'^^ , y = l/VlTI • (33) 

(2) Fitting functions to the DE density 

The dark energy density can be written as a truncated Taylor series polynomial 
in X = 1 + z, pde = Ai + A2X + Asx^. This leads to the following ansatz for 
the Hubble parameter 

H{x) = Ho[fl^x^ + Ai+A2X + A3X^]^ , (34) 
which, when substituted in the expression for the luminosity distance (9), yields 
c dx ^3^^ 



1 + z Ho Ji V^mX^ + ^1 + A2X + A-iX^ ' 



Reconstructing Dark Energy 11 



This fitting function gives exact results for the cosmological constant w = —1 
{A2 = A3 = 0) as well as for quiessence with w = —2/3 {Ai = A3 = 0) 
and w = — 1/3 {Ai = A2 =0). The presence of the term il,„x'^ in (34) ensures 
that the ansatz correctly reproduces the matter dominated epoch at early times 
{z ^ 1). For quintessence models as well as the Chaplygin gas, the luminosity 
distance Dl{z) in (35) can be determined to an accuracy of better than 1% if 
n^o > 0.2 38. 
(3) Fitting functions to the equation of state 

In this approach one assumes that the DE equation of state w{z) is an unknown 
variable whose behaviour is 'guessed' by means of a suitable fitting function 
w{z,ai). Since 



Six = (1 - ^^mo)exp|3y j-^p^ — -dz^ , 



(36) 



the values of a; can be determined by comparing the luminosity distance (9) or 
D A{z)^r{z), against observations. 

Several possible fits for w{z) have been suggested in the literature. Perhaps the 
simplest is the Taylor expansion 

N 

w{z) = ^Wi^;' , (37) 

i=0 

which, for A'' = 1, gives results significantly better than the Taylor expan- 
sion for the luminosity distance (30). The simple two parameter representation 
w(z) = wq ^ w\z is however of limited utility since it is only valid for 2; <C 1. A 
considerably more versatile four parameter ansatz has been suggested in . 
A popular two parameter fit was suggested by Chevallier and Polarski and 
by Linder 

w{a) = + w\{\ - a) = -\- w\- , (38) 

1 + -2 

the associated luminosity distance can be obtained by substitution into (36) 
and (9). A more general form for this fit is 

w{a) = + (ttp - a)Wa , (39) 

where Op is the value of the scale factor at the 'sweet spot' where the equation 
of state w(a) is most tightly constrained. The value of usually depends upon 
the data set being used 87,88,89,90 

As mentioned earlier, for an ansatz to be regarded as being successful it should 
embrace within its fold the behaviour of a reasonably wide class of DE models. 
Most of the ansatz's discussed above successfully accommodate DE whose equation 
of state evolves moderately with redshift. It is quite clear that these simple fits 
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Gold= 1 56 SNe _ 




Fig. 1. An early reconstruction of the supernova data (left panel) using the parametric fit 
(34) shows an evolving equation of state to be marginally preferred over the cosmological 
constant The same data were independently analyzed by means of a non-parametric 
ansatz (right panel) . It is encouraging that both reconstructions appear to give similar 
results for Q < z <1 where most of the data points lie. The crossing of the so-called 'phan- 
tom divide' at ui = — 1 led to much theoretical interest and considerable model building 
activity. More recent SNe results, together with constraints from other measurements such 
as the CMB, LSS and Baryon acoustic oscillations imply less pronounced evolution of w{z) 
with redshift as shown in figure 2, see also 49,50,51,79,82,83,91,92,93,94 r^j^^ j^^^ 

from Alam, Sahni and Starobinsky while the right panel is from Huterer and Coorey 
83 



cannot be used to rule out (using observations) models with rapidly evolving 'w{z). 
One reason for this is that applying an ansatz such as (34) or (38) to SNe data is 
equivalent to smoothing the evolution of the Universe over a redshift interval Az oc 
1/N where A'^ is the number of free parameters in the ansatz (the concrete form of 
an ansatz defines the form of a smoothing redshift function). Clearly, an implicit 
smoothing such as this will cause rapid transitions in dark energy to disappear not 
because of disagreement with data but simply because of the manner in which the 
data have been 'massaged' 39, 40^ (^p^^ example of the disastrous results of applying 
the prior w — constant to models with an evolving EOS is discussed in section 5.) 
Models with a fast phase transition in dark energy (|^| ^ 1 over a narrow range 
of redshift 5z <C 1) have been discussed in ^8,59_ accommodate such models the 
following ansatz was suggested ^^'^'^ 

"^^^-"'+ l-^exp(^) ' ^''^ 
where Wi is the initial equation of state at high redshifts, zt is a transition redshift 



Reconstructing Dark Energy 13 



at which the equation of state falls to ■w{zt) = {wi + Wf)/2 and A describes the rate 
of change of w{z). Substitution into (36) and (9) gives the luminosity distance. 

Finally, models with oscillating DE have also been discussed in the literature 
60,49,52,54,94_ -j-q accommodate an oscillating EOS the following ansatz has been 
suggested 

w(loga) = wo+wi cos [Aloga/ad (41) 

where wq, wi,A, Oc are free parameters whose values must be obtained by fitting to 

observations. It is encouraging to note that, in spit(^ of some ambiguity in the form 
of the different fits, when applied to the same supernova data set most of them 
give consistent results in the range 0.1 < z < 1 (where there is sufficient data) 

39,40,49,51,83 

One should also note that although increasing the number of parameters usually 
increases the accuracy of reconstruction of the 'best fit', this is often accompanied 
by severe degeneracies which limit the utility of introducing a large number of 
free parameters. In addition, extra free parameters are quite severely penalized 
by information criteria such the Akaike information criterion and the Bayesian 
information criterion see also ^8,63,64 However, there exists a subtlety when 
applying these information criteria to fitting functions of H{z), w{z) etc. as opposed 
to concrete theoretical models of DE. In the case of the former, one should keep in 
mind the possibility that a fit with a larger number of parameters might in fact be 
describing the behaviour of a fundamental DE model (not yet known) containing a 
smaller number of truly free parameters. 

The reverse is also triie, a primitive fit (hence not penalizablc) may. if the in- 
formation criteria are applied, detract our attention from a more complicated but 
also more fundamental explanation of a phenomenon. Let us illustrate this with an 
example. The theory of gravitational instability informs us that the linear density 
contrast 6 and the peculiar velocity field v are related as ^^'^^ 

^ = -^V.v (42) 

where the function / = (ilog(5/(iloga can be approximated as / ^ ^rif- An exact 
calculation however reveals / to be a more complicated function involving elliptical 
integrals, etc. Therefore a naive application of information criteria to (42) would 
needlessly penalize the latter, which is correct but complicated, in favour of the 
simpler / ~ fl^. 

Another popular approach the so-called Principal Component analysis is based 
on expanding u){z) in terms of a basis of orthogonal functions. The form of these 
functions depends upon the kind of data used and its constraining capabilities 
67,68,83 ^j^j ^Yiis section by noting that the reconstruction approach has also 
been extended to scalar-tensor gravity ^, string inspired cosmology with the Gauss- 
Bonnet term coupled to a scalar field and dark energy with non-canonical 
kinetic energy terms Other applications of cosmological reconstruction may be 
found in 
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Fig. 2. The equation of state of dark energy w{z) reconstructed using the WMAP 3 year 
data + 157 "gold" SNIa data + SDSS. Median (central line), 68%(inner, dark grey) and 
95%(outer, light grey) intervals. The two parameter fit (38) has been used in this exercise. 
From Zhao et al. . 

4. Non-parametric reconstruction 

Non-parametric smoothing usually involves directly smoothing either D^, or some 
other quantity appropriately binned in redshift space with some characteristic 
smoothing scale. Different ways of implementing this approach have been discussed 
73.67,74.75,76.78,79.80.81,82,83,84 All reconstruction methods must deal with the fact 
that the data sample is usually sparse and the coverage of redshift space uneven. 
Consider for instance the comoving coordinate distance to a single supernova within 
a larger sample 

n = r{z,) + n.i . (43) 

Here, rii is the noise term {rii) = 0, {niUj) = crfSij, and r(z) is related to the distance 
modulus /io of the supernova by r(z)/lMpc = 10^0/^-^/2997.9(1 + z). Direct 
differentiation of to give H~^{zi) will obviously result in a very noisy quantity. 
One way to tackle this is to fit ri[z) piece- wise using a set of basis functions. 
Following this approach, Daly and Djorgovsky ^^'^^ reconstructed the deceleration 
parameter from the dimensionless coordinate distance y = Hqv by: (i) fitting y{z) 
locally within a small redshift bin by means of a second order polynomial requiring 
that at least 10 data points lie within each bin; (ii) determining the first and second 
derivatives of y{z) from the fit coefficients (for each bin) and reconstructing the 
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deceleration parameter from y{z) by means of the relation 

^ = 1 + (1 + .)^. (44) 

Applying their method to a set of radio galaxies in addition to type la SNe, Daly 
and Djorgovsky 75.76,77 found that although the cosmological constant provided a 
very good fit to the data, modest evolution in the dark energy density was also 
perfectly acceptable. The transition from deceleration to acceleration occured at 
z > 0.3, with a best fit value of ^; = 0.42. These results are in broad agreement with 
those obtained using parametric approaches 39.40, 49,51 

A different approach to non-parametric reconstruction is discussed by Shafieloo 
et al.^^ who generalize a smoothing ansatz widely used in the analysis of large scale 
structure. According to this method, a smoothed quantity D'^(x) is constructed 
from a fluctuating 'raw' quantity -D(x') using a low pass filter F having a smoothing 
scale A 

D^(x, A) = y" -D(x')F(|x - x'l; A) dx' . (45) 

In large scale structin-c studies, D is the density field ^^'^^ whereas for cosmologi- 
cal reconstruction D could be either of Dl{z), DA{z), r{z). Commonly used filters 
include: (i) the 'top-hat' filter, which has a sharp cutoff 



^TH 



a e 1 



where 9 is the Heaviside step function (©(2:) = for 2; < 0, 6(2;) = 1 for 2; > 0) 
and (ii) the Gaussian filter 



Fg (X exp ' J - (46) 

When applied to SNAP-quality data using a Gaussian filter, this method recon- 
structs the Hubble parameter to an accuracy of <, 2% within the redshift interval 
< 2 < 1. The look-back time 

/■^ dz' 

is reconstructed to an even better accuracy of <, 0.2% at z ~ 1.7. 

Despite the considerable success of these approaches, neither they, nor the para- 
metric methods discussed earlier provide independent measurements of either H{z) 
or q{z) within a given redshift range. One approach towards an uncorrelated de- 
termination of cosmological quantities is to bin the data in redshift bins of size A 
with a weight (or window) function assigned to each bin Each window function 
vanishes outside its redshift bin and the two supernovae falling on either side of a 
bin boimdary arc discarded. This method ensures that, since two adjacent window 
functions have no common supernovae, the measurement of quantities within these 
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Fig. 3. Dark energy density px{z) reconstructed using SN la data combined with 
the WMAP 3 year data, the SDSS baryon acoustic oscillation data, and the 2dF linear 
growth data for density perturbations. The 68% (shaded) and 95% confidence contours are 
shown. Beyond 2cut = l-4 (upper panel) and f.Of (lower panel), px{^) is parametrized by 
a power law (f + z)". The horizontal dashed line shows the unevolving density associated 
with a cosmological constant. From Wang and Mukherjee 

bins will have uncorrelated error bars. This method is quite successful and com- 
plements other reconstruction approaches in many respects. One weakness, which 
practitioners of this approach have to guard against is that, since some information 
at the bin boundary is lost, discontinuities in the fit values at bin edges may be 
present These, in turn, could lead to unphysically large derivatives of derived 
quantities such as H{z) and u;(z). Results obtained by applying a variant of this 
method to recent CMB+SNe+LSS data are shown in figure 3 and appear to agree 
with the results of the parametric approach in figure 2. Other promising methods 
for obtaining uncorrelated estimates of cosmological evolution are discussed in 

In all of the above approaches, it is important to choose the value of the smooth- 
ing scale A (or bin size) optimally, so that (i) a sufficient number of data points is 
accommodated within each bin to reduce the effects of shot noise, (ii) the bin size 
is not too large to cause excessive smoothing. The following useful formula giving 
the relative error bars on H{z) comes to our aid 
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where N is the total number of supernovae (assuming an approximately uniform 
distribution) and a is the noise of the data. Clearly, a reduction in the value of A 

by 3 should be compensated for by having 27 many times more supernovae if one 
wants to keep ^ unchanged. The situation is far worse for w{z) which has a A"^/^ 
scaling, since the data are being differentiated twice to get 'w{z). 

5. Obstacles to cosmological reconstruction 

With respect to cosmological reconstruction, considerable caution must be exercised 
when setting priors on the values of cosmological parameters. This is true both for 
the density parameter whose value is currently known to about 15% accuracy 
and the DE equation of state. Maor et al^^ gave a particularly insightful example 
of the dangers of incorrect reconstruction by assuming a fiducial DE model with 
= 0.3 and an evolving EOS wq{z) = —0.7+ 0.8z. In their reconstruction, the 
results of which are shown in figure 4, it was assumed that wq was a constant and, 
additionally, that wq > — 1. These two (incorrect) priors led to a gross underesti- 
mation of Wq and an overestimation of il„, as shown in figure 4. 

The above example shows that the EOS of dark energy can be badly recon- 
structed if we assume a misleading prior for w. Since w{z) depends upon the value 
of ri„i through (7), it follows that an incorrect assumption about the value of the 
matter density has the potential to affect cosmological reconstruction quite signifi- 
cantly. This is demonstrated in figure 5 in which the DE EOS is reconstructed for 
a fiducial ACDM model with w{z) = — 1 and fim = 0.3. Cosmological reconstruc- 
tion is based on the ansatz (34) in which the incorrect value flm = 0.2 is assumed. 
This results in the reconstructed EOS evolving with redshift when, in fact, no such 
evolution is present in the fiducial model ! 

6. Looking beyond the equation of state 
6.1. The w-probe 

As we have seen, the EOS has certain blemishes - it is obtained from Dl after 
differentiating twice and it is sensitive to the precise value of fimo- Both these 
features hamper its reconstruction which has led several authors to propose the 
Hubble parameter or the DE density as being more appropriate for cosmological 
reconstruction 26,75,76,47,79_ Carrying these ideas further, the possibility of extracting 
information about the equation of state from the reconstructed Hubble parameter 
by constructing a weighted average of the equation of state was explored in This 
quantity, dubbed the u>-probe, gleans information about the equation of state from 
the first derivative of the luminosity distance. Therefore, it is less noisy and better 
determined than w{z). 

The weighted average of the equation of state is defined as 




(49) 
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Fig. 4. Pitfalls in cosmological reconstruction 1. This figure from Maor et al. shows 
the results of a reconstruction exercise performed by assuming: (i) wq = constant as a 
prior, which gives the larger lower contour with wq < — 1; (ii) the additional constraint 
WQ > — 1, results in the smaller upper contour with wq = — 1 as the best fit. Both (i) and 
(ii) give confidence contours and best fit values of wq and Qm which are widely off the 
mark since they differ from the fiducial Quintessence model which has wq{z) = —0.7+0.80 
and fim = 0.3. 



An important feature of the w-probe {w) is that it can be expressed in terms of 
the difference in dark energy density over a given redshift range 



l + w{Zi,Z2) = 



1 S InppE 
3 6 ln(l + z) 



lln 
3 



H^{Z2)-n0m{l+Z2)'' 



'In 



l + zi 

1 + Z2 



(50) 



Here 5 denotes the total change of a variable between integration limits and pde = 

Pde/ pac {pqc = SHq/SttG). Thus w is easy to determine if the Hubble parameter 
has been accurately reconstructed. 

An important property of iu is that it is less sensitive to uncertainties in the 
value of OmO than w{z). This is demonstrated in figure 6 where the results for Hj are 
shown after marginalizing over the matter density. Remarkably, the value of iij for 
the fiducial ACDM model remains close to —1, while w for the evolving DE model 
shows a clear signature of evolution. Thus, small uncertainties in the value of the 
matter density do not appear to adversely affect the accuracy of the reconstructed 
w-probe. Furthermore, several excellent methods for determining pde and H{z) 
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Fig. 5. Pitfalls in cosmological reconstruction 2. The reconstructed equation of state ■w{z) 
is shown for 1000 realizations of an Q„i — 0.3, w — —1, ACDM model assuming SNAP 
quality data. An incorrect value for the matter density, fimo = 0.2, is assumed in the 
reconstruction exercise which uses the polynomial ansatz (34) . The dashed line represents 
the fiducial ACDM model with w = —1 while the solid lines show the mean value of the 
(incorrectly) reconstructed w(z) and la confidence levels around the mean. Note that the 
reconstructed EOS excludes the fiducial ACDM model to a high degree of confidence. From 
Shafleloo et al. 

have been suggested in the Hterature 75,76,26,39,47,79^ ^^^y of which could be used to 
determine w using (50). 

6.2. A geometrical diagnostic of dark energy 

Cosmological observations made during the past two decades have brought about 
both qualitative and quantitative changes in our perception of the Universe. Obser- 
vational support has come for many theories (such as inflation) which earlier were 
regarded as being at the level of hypothesis. On the other hand, fundamentally new 
properties of our Universe, such its current state of acceleration, have also been un- 
ravelled. The remarkable qualitative similarity between current acceleration (fueled 
by DE) and inflationary acceleration (supported by what may be called primordial 
dark energy) may help us understand in which direction cosmology may be headed. 
The discovery by COBE of an approximately scale-invariant spectrum of primor- 
dial fluctuations (n^ — 1 ~ 0) was regarded by many as providing tacit support for 
the inflationary scenario. However, virtually all models of inflation predict small 
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Fig. 6. The w-probe is reconstructed for the ACDM model with w = ~1 (left) and an 
evolving DE model with w = — l/(l + z) (right). 1000 realizations of SNAP- like data have 
been used. The thick dashed line in both panels indicates the exact value of w for the 
fiducial model, the dark grey boxes in each panel indicate the Icr confidence levels on w 
reconstructed for the two models after marginalising over = 0.3 ± 0.07. The results 
of both panels show that the w-probe is able to accurately determine DE properties even 
if the matter density is not perfectly known. Prom Shaiieloo et al. . 

departures from the scale invariance Ug — 1 — where e is a model dependent 
quantity. Strict scale invariance [ris = 1) arises for a specific 2-parameter family of 
infiaton potentials only which reduces to V{<t)) oc in the slow-roll approx- 
imation. In this context, the recent 3-year data release of the WMAP experiment 
suggests Us — 1 ~ —0.04, which is in excellent agreement with predictions made by 
the simplest inflationary models such as chaotic inflation with V{(j)) cx 0^. 

The present situation concerning DE resembles, in some respects, the status 
of inflation just after the release of the COBE data in 1992. The cosmological 
constant is in excellent agreement with data and departures (if any) from 1 + = 
are believed to be quite small. Nevertheless, since at present no fundamental 
theory predicting the (small) value of the A-term exists, it is of utmost importance 
that departures from 1 + w = Q he probed to accuracies of at least 1% by future 
generations of dark energy experiments. Since virtually all DE models (other than 
A) have either wq ^ —1 ov wq ^ the need for probing both the equation of state 
as well as its first derivative become crucial if any deep insight is to be gained into 
the nature of DE. 

It is useful to recall that several DE models which agree well with the cur- 
rent data arise because of modifications to the gravitational sector of the theory 
7,8,io,96,97_ Pqj. ^jjggg geometrical DE models, the EOS no longer plays the role of a 
fundamental physical quantity and it would be very useful if we could supplement 
it with a diagnostic which could unambiguously probe the properties of all classes 
of DE models. The w-probe discussed earlier, provides one such method, since it 
is based on the expansion history H = a/a. Since the expansion factor a{t) is an 
essential feature of all metric theories of gravity, it is worthwhile trying to explore 
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the properties of DE by considering the generic form 



a{t) = a{to) + d|o(i - to) + -^{t - tof + -^(t - + ... . 



(51) 



Cosmic acceleration appears to be a fairly recent phenomenon 98,99,39,40^ gg ^g^j^ 
confine our attention to small values of |i — to| in (51). The second and third terms in 
the RHS of (51) have played a crucial role in the development of cosmology. Indeed, 
for a long time the large errors in H = a /a and q = —d/aH^ impaired a precise 
picture of cosmic expansion and only recently has the sign of q boon determined to 
sufficient accuracy for us to make the statement that the universe is accelerating. ^ 
Keeping in mind the significant progress expected in observational cosmology over 
the next decade, we feel the time to be ripe to supplement H and q with r = d /aH^, 
which is the next logical step in the hierarchy of cosmological parameters. It is 
remarkable that r = 1 for the spatially flat ACDM model which is the 'standard' 
cosmological model at present, as well as for the CDM model. Supplementing r (the 
first statefinder) with the second statefinder s = (r — l)/3(g — 1/2) permits us to 
break this degeneracy and to characterize different DE models in a very informative 
manner. ^ 

The statefinder pair {r, s} is a 'geometrical' diagnostic since it depends upon a{t) 
and hence upon the space-time geometry. An important property of the statefinder 
is that spatially flat LCDM corresponds to the fixed point 



So the basic question of whether DE is the A-term or 'something else' can be 
rephrased into whether the equality (52) is satisfied for DE. Indeed, the departure 
of a DE model from {r,s} = {1,0} provides a good way of establishing the 'dis- 
tance' of this model from LCDM 26,38_ -j-j^g jg^^^ ^j^^^ different classes of DE models 
show distinctly different behaviour when plotted in the {r, s}-plane adds to the 
practical utility of this diagnostic. See figure 7 and ^s.ioi £qj. gome applications of 
the statefinder. 

Finally, the pair {r, s} can also be expressed in terms of {w, w} quite simply 



{r,s} 



= {1,0}. 



(52) 



LCDM 



, 3^ ti; 

1 w 




(53) 



(54) 



It is instructive to note that the second statefinder can be rewritten as 



38 



s = 



{P + P) P 
P P 



(55) 



'The geometrical relation R/6H^ = r2total — <7 links the Ricci scalar R to q. Therefore, an accel- 
erating (q < 0) spatially flat universe corresponds to R/QH^ > 1. 
sThe first statefinder r was also dubbed the cosmic jerk in 1*^°. 
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Fig. 7. The time evolution of the statefinder pair {r, s} for quintessence models and the Chaplygin 
gas. Solid lines to the right of LCDM represent the inverse power-law potentials V = Vo/<l>", 
while those to the left correspond to the Chaplygin gas. Dot-dashed lines represent DE with a 

constant equation of state w. Tracker models tend to approach the LCDM fixed point (r = 1, s = 0) 
from the right at t — > oo, whereas the Chaplygin gas approaches LCDM from the left. For Chaplygin 
gas K is the ratio between matter density and the density of the Chaplygin gas at early times 
The dashed curve in the lower right is the envelope of all quintessence models, while the 
dashed curve in the upper left is the envelope of Chaplygin gas models (the latter is described 
by K = Hm/l — ^m)- The region outside the dashed curves is forbidden for both classes of dark 
energy models. The ability of the statefinder to differentiate between dark energy models is clearly 
demonstrated. From Alam et al. 

where p = Xla^'o total pressure including that of DE. From (55) we find that 

s is very sensitive to the epoch when the total pressure in the Universe vanishes. 
For ACDM, this takes place at the rcdshift Zp ~ 10, when the positive radiation 
pressure cancels the negative pressure in A. Consequently, s ^ oo as z ^ Zp. In 
general, the redshift Zp at which the total pressure in the Universe becomes zero is 
quite sensitive to the nature of DE. 

Finally, returning to our analogy between an inflaton and dark energy, an in- 
teresting correspondence exists between the {q, r} pair describing DE and the slow 
roll parameters e = —H/H^, rj = —H/2HH, describing inflation 



q + l = e 

r - 1 = e (2?? - 3) 



(56) 
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Level Geometrical Parameter Physical Parameter 

1 H{z)=l^ Pm{z) = POm{l + zf, 

PDE = 1^ — Pm 



_1 I rflogg 
^ d log (1+0) 



ACDM 



3H-' 



2 ' "'^"^ T+V 
O _ 8-kGV o — SttGT 



r-l 
3(9-1/2) 



{r, s} 



{1.0} 



ACDM 



We end this section by noting that our current description of cosmology relies 

on parameters which arc cither geom,etrical or physical in nature. In the tabic above, 
we have divided cosmological parameters according to their level which is related 
to the number of differentiations of the expansion factor needed to construct that 
parameter. Note that all relative energy densities in the tabic Vla{z) are defined 
using the present value G = Gq of the Newton gravitational constant in agreement 
with the definition (1). Each higher level requires an additional differentiation of 
observational data {Dl,Da, ■■■) and, therefore, demands a higher level of accuracy 
for the latter. For current data, we broadly have about 10% accuracy for quantities 
belonging to the first level, 50% for the second level, while determination of third 
level parameters lies in the future. By smoothing, we lower the level of a parameter 
by unity. 

Relations between geometrical and physical parameters are the following: 



^T{z) = ^^-^nrno{l + zf ; 

nn{z) = -i^r-3q-A+y^nmo{l + zfj . 
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7. Summary and discussion 

The nature of dark energy is clearly one of the outstanding physical and cosmolog- 
ical puzzles of this century. Although many distinct theoretical models have been 
advanced to explain the cosmic acceleration, an alternative approach to study and 
understand DE is to determine DE properties from observational data in a model 
independent fashion, thus, reconstructing DE from observations. In this review, we 
have attempted to briefly summarize some important approaches of cosmological 
reconstruction. These approaches are briefly classifiable into parametric and non- 
parametric methods. Both methods appear to give consistent results when applied 
to current data. Future directions for model independent reconstruction of dark 
energy have also been briefly discussed. In the first approximation, within current 
observational errors, the reconstructed DE behaves like a cosmological constant. But 
future, much more precise data, may well show some deviations from this behaviour. 
However small, these could have a profound influence on the whole of physics. 
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